(f) production at RHIC: characterization of co-existence phase 
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We extract the effective degrees of freedom that characterize the co-existing phase of quark gluon 
plasma and hadrons. Experimental data on <f> at mid-rapidity is used to set a lower bound to the 
critical temperature of quark hadron phase transition. The production and evolution of strangeness 
have been studied by using Boltzmann equation. The results have been contrasted with the experi- 
mental data obtained by STAR collaboration at RHIC for Au + Au collisions at ^/sjvjv = 200 GeV. 
Our study reveals that the 4> mesons freeze out at a temperature ~ 160 MeV, a value close to the 
transition temperature for quark-hadron phase transition. 
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The goal of the ultra-relativistic heavy ion col- 
lisions is to create and study quark gluon plasma 
(QGP). The pre-requisite for the formation of QGP 
is to create a hot and dense system of nuclear matter. 
Experimental results from the Relativistic Heavy- 
Ion Collider (RHIC) have indicated the formation of 
such a state of matter in the early stage pj . Among 
others, one of the promising signals of QGP forma- 
tion is the enhanced production of strange and anti- 
strange particles Q. As a result, the production 
of particles with hidden (e.g. 4>) as well as open 
strangeness have become a field of intense theoreti- 
cal [iHHEi and experimental H0HE3 activities - 



In the present work we study the production of 
(j> meson at RHIC. The <fi meson is composed of a 
strange (s) and anti-strange (s) quark and its inter- 
action with nuclear matter is suppressed according 
to Okubo-Zweig-Izuka (OZI) rule. In QGP s and s 
quarks are produced mainly by gluon fusion and an- 
nihilation of light (u and d) quarks and anti-quarks. 
These ss will form (f> through the hadronization pro- 
cess. This process is not OZI suppressed. There- 
fore, we expect excess 4> if QGP is formed in the 
initial state of heavy ion collisions as compared to 
hadronic initial state. The production of <j) mesons 
from sources other than plasma is expected to be 
small, so they will not overshadow the plasma signal. 
The other advantage of d> is that after its produc- 
tion during the hadronization of the QGP, it suffers 
less re-scattering in the hadronic matter, thereby re- 
taining information of the thermodynamic state of 
the matter during its hadronization stage. There- 
fore, study of cj) spectra will be very useful to de- 
termine the transition temperature (T c ) for quark- 
hadron phase transition. However, it may be men- 
tioned that the kaons which come from decay of tj> 
mesons will undergo re-scattering, and the </>'s exper- 
imentally reconstructed from these kaons may carry 
some re-scattering effect. Still reliable information 
on T c can be extracted by measuring 4> productions 
in high energy heavy ion collisions [llj. 



Hence we shall show the usefulness of <j> as a probe 
for co-existence phase by studying its production at 
y/SNN = 200 GeV First we demonstrate from 
the experimental data that <fi meson freezes out at 
a temperature close to T c , the value of which is ob- 
tained from lattice QCD simulations 12]. We argue, 
using the yield of 4> mesons at mid-rapidity (-g-jj^) 
that the enhanced production is possible if the sys- 
tem has passed through a mixed phase of quarks, 
gluons and hadrons and the temperature during the 
hadronization is > 160 MeV. Finally, attempt will 
be made to characterize the mixed phase by putting 
bounds on the the effective statistical degeneracy. 

Relativistic hydrodynamics with (3+1) dimen- 
sional expansion |13| and boost invariance along the 
longitudinal direction 0] has been used to study 
the evolution of the matter formed in heavy ion col- 
lisions. We take the equation of state from lattice 
QCD simulations 01 • The initial velocity and en- 
ergy density profiles are taken as v(r,Ti) — and 
e(r, T}) = eo/(l + exp(r — R)/S) respectively. Here Tj 
is the thermalization time, R is nuclear radius. Re- 
sults will be shown for two sets of values of r, = 0.15 
and 0.2 fm/c. Such values of r, 0, 0] have re- 
cently been used to reproduce the thermal photon 
spectra measured by PHENIX collaboration 0] in 
Au + Au collisions at ^/s~n~n = 200 GeV. The cor- 
responding values for the initial temperatures are 
590 and 480 MeV respectively. Results of trans- 
verse momentum spectra of <p mesons are shown 
in fig. ^ It appears that the results for T, = 590 
MeV describes the data quite well for the entire 
range of transverse momentum (j>t) under consid- 
eration. However, for px more than 2 GeV (or 

rriT — ~ 1.2 GeV, where rax — \Jp T + ^ s 
the transverse mass and is the mass of the 4> me- 
son) the application of hydrodynamics may be ques- 
tionable and contributions from perturbativc QCD 
may become important. In fact, at low px domain 
the contribution of thermal ss for 4> formation is 
dominant 18] . In that case, the description for 



2 




FIG. 1: Transverse mass distribution of cj> mesons for 0- 
5% central Au+Au collisions at ^/snn = 200 GeV. Com- 
pared to hydrodynamical calculations with and without 
radial flow and for two different initial conditions. The 
equation of sate used is obtained from lattice QCD re- 
sults. 



FIG. 2: The time evolution of the ratio of non- 
equilibrium to equilibrium number density of strange 
quarks/hadrons for transition temperature T c — 150 
MeV. The variation of temperature with time is also 
shown. 



Ti = 480 MeV is better suited. From the inverse 
slope (for 0.45 < p T (GeV)< 2) of the data we get 
the effective temperature as 257 ± 8 MeV. We find 
that the average radial velocity < v r >~ 0.43 at 
the freeze-out surface. Therefore, using the relation 
T e ff — T true + (l/2)m0 < v r > 2 , we get the 'true' 
freeze-out temperature T true ~ 162 ± 8 MeV (error 
bars correspond to the experimental error bars @)- 
It should be noted that the value of T true is close 
to T c predicted by lattice QCD simulation. The 
value of the hadronization temperature in statistical 
hadronization model 0] was found as T c = 155 ± 8 
MeV (see also |20|). The similar value of hadroniza- 
tion temperature was obtained |2l| from the analy- 
sis of experimental data on heavier hadrons. 

Now we turn to the yield of </> mesons at mid rapid- 
ity. The experimental value of for 0-5% central 
collisions is 7.7 ± 0.30 |8(. It is more convenient to 
deal with the ratio, -jjf-/^-, where dN/dy is the 
total particle multiplicity at mid-rapidity. This ra- 
tio has the advantage in the sense that some of the 
uncertainties (e.g. volume)are removed. From ther- 
mal model calculations it can be shown that at the 
critical temperature T c , 

where, 771^=1.02 GeV, and <? e // is the effective 
statistical degeneracy. Eq.^ indicates that the ratio 



~~d^~/(i!y depends on T c and the statistical degener- 
acy where the </> mesons freeze out. Taking T c and 
g e ff from lattice QCD simulations ^3] as 170 MeV 
and 18 respectively, we get = 4.77 x 10~ 3 

from Eq. ^ The experimental value of -j^/t^t is 
7.7xl0~ 3 . This indicates that there is a clear over 
production of <p mesons. The data can be expressed 
in terms of the calculated value of dN$ /dy \ eq by a 
simple relation: 

_ dN^/dyl 

expt /n\ 

70 = dN/dy\ eq [Z) 

dN^/ dy\ eq is the value of <fi meson multiplicity ex- 
pected from thermal model at T c and 7^ is the en- 
hancement factor (a value of 7^ = 1 (< 1) indicates 
equilibrium (under production). The experimental 
data requires the value of 7^ to be ~ 1.3. This in- 
dicates that there seems to be a over saturation of 
strangeness production at the temperature where <f> 
freezes out. 

Now we ask, at what stage of the evolution of the 
system formed in heavy-ion collisions gives the re- 
quired over production of strangeness? We adopt the 
momentum integrated Boltzmann equation to study 
the evolution of the strangeness density. This has 
also been used extensively to investigate the density 
evolution in the early universe |22j | . The initial state 
is considered as a thermalized state of gluons and 
light quarks with initial temperature, ~ 480 MeV 
and thermalization time, Tj ~ 0.2 fm/c. The reac- 
tions considered for the strange quark productions 
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FIG. 3: Same as FigHfor T c = 170 MeV. 
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FIG. 4: Same as FigEJfor T c = 195 MeV. 



in QGP are: qq — > ss and (75 — > ss. For the produc- 
tion of strange hadrons the following reactions are 
considered 123: tt + + n~ — > K + + K~ , 7r + + tt° — > 



Jf + + K° and 7T° + tt° 
meson sector and 7r + +n 



+ in the meson- 

p—> Y.+ +K+, 7r~ T,~+K + and 7r° + n 

— » £~ + for the meson-baryon sector. 

The evolution of s in QGP phase with proper time 
r is given by 



dr s 
~d^~ 



eq 



(3) 



Similar equation can be written down for s quarks. 
In the present case r s = r s . The evolution of K + (us 



mesons in the mixed phase is governed by the fol- 
lowing equation: 



dr K + 
dr 



Wq (1 - r K +r K - 

— eq — (1 ~ r K +r A ) 

eq 
n K+ 



(1 - r K +r^) 



ldf ( r- s n 



eq 



f dt \ 



(4) 



In the above equation the last term stands for 
the hadronization of s quarks to K + 0, HH- ol 
is a parameter which indicates the fraction of s 
(s) quarks hadronizing to K~_(K + ), the value of 
a = 0.5 if we consider K + , K a , K~ and K° for- 
mation in the mixed phase, denotes the ratio of 
non-equilibrium (rii) to equilibrium (n^ q ) density of 
the particle i. Similar equation can be written for 
rfc in the hadronic phase. The evolution of other 
particles e.g. K~ can be treated similarly with ap- 
propriate cross sections [23^ as inputs and solving the 
coupled set of equations. The initial value of s and 
s quarks are taken close to their equilibrium value. 
However, a small change in the initial values of r s do 
not change the final results drastically. Even with 
lower initial values of s and s the system reaches 
equilibrium very fast due to their production in the 
high temperature heat bath. In Eq. 0|/ and (1 — /) 
are the fractions of hadrons and QGP respectively 
in the mixed phase. Ri(T) is the production rate of 
particle i at temperature T. The cooling of the heat 
bath is governed by the following equation [26j: 



dT 
dr~ 



b(r s 



rs 



a(a + b(r s + r s )) 



(5) 



where a = (1 + c 2 s )/c 2 s , a = 8tt 2 / '(45c;?) and b = 
77T 2 np'/(120Cs), c 2 is the velocity of sound, np is the 
number of flavours and r s — dr s /dr. It may be noted 
that the Bjorken's cooling law can be recovered from 
Eq. El when the second term vanishes. The value of 
net baryon number at the central rapidity has been 
assumed to be zero. 

The time evolutions of the ratios of non- 
equilibrium to equilibrium number density for 
T c = 150, 170 and 195 MeV are shown in Fig. 03 
and0]respectively. It should be mentioned here that 
the corresponding g e f t for the above values of tem- 
perature are taken from lattice QCD [l2|. We ob- 
serve a clear over saturation in number density of 
kaons at the end of the mixed phase (Fig. |3k0J. 
For T c = 170 - 195 MeV the value of r K is 1.4 - 1.5 
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at the end of the mixed phase, before it reaches the 
equilibrium value of 1 in the hadronic phase. It may 
be noted that there is a small difference between 
and r K - in the hadronic phase because of their dif- 
ferent production cross sections. For the value of T c 
= 150 MeV, we observe that ~ 1, for a = 0.5. 
For smaller values of a, r K + will be even smaller. 

Considering that we need a 7^ > 1 to explain the 
experimental data (as explained before) our calcula- 
tions set a lower bound on the value of T c , it must 
be > 160 MeV. For a hadronic initial state with 
temperature of 250 MeV the value of Tk remains 
< 1. This indicates that the enhanced production 
of (j> mesons in heavy-ion collisions at ^snn = 200 
GeV can be explained if the system passes through 
a mixed phase during it's evolution with a T c > 160 
MeV. At this value of T c , the enhancement factor, 
70 ~ r K+ x r K- ~ Ik ~ 1>24. Considering the fact 
that the experimental data indicate 7^ ~ 1.4 ± 0.1 
we argue that T c > 160 MeV. 

The life time of the mixed phase crucially de- 
pends on the effective statistical degeneracy of QGP 
and hadronic phases. g e tf ~ 9,26 and 31 for T = 
150, 170 and 195 MeV respectively. Larger difference 
in the statistical degeneracy gives rise to larger dis- 
continuity in the entropy density which consequently 
generates more latent heat. In the mixed phase the 
cooling due to expansion is compensated by the lib- 
eration of latent heat. Hence if larger latent is gen- 
erated then the system will spend more time in the 
mixed phase i.e. the life time of the mixed phase 
will be larger. This is clearly reflected in the results 
shown in Fig. 121 II 

We characterize the mixed phase in terms of ef- 
fective statistical degeneracy g e ff- This is justified 
because, (i) <j> mesons freeze out at T c and (ii) its 
production can only be explained by the existence 
of mixed phase. For fixed g e ff we calculate -js^/^p 
from Eq. ^ This ratio is then multiplied by 7^ and 
compared with the corresponding value from exper- 
imental data (lies within the two thick horizontal 
lines in Fig. EJ. The results for g e ff — 18, 24, and 
30 are shown in Fig. EJ We find that g e ff ~ 18 may 
be considered as a lower bound. 



In summary, we have studied the production of <j> 
mesons in Au+Au collisions at ^Jsmn — 200 GeV. 
The transverse momentum spectra is well explained 
by hydrodynamical model with the equation of state 
taken from lattice QCD. It is observed that </> mesons 
freezes out around the transition temperature pre- 
dicted by lattice QCD. The over production of <p 
meson at mid-rapidity can only be explained if QGP 
is formed in the initial state and passes through a 
co-existence phase of QGP and hadrons, i.e. over 
saturation in strangeness indicates the formation of 
co-existing phase of QGP and hadrons (indications 
of which has already been pointed out at SPS en- 
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FIG. 5: The variation of (dN r/> /dy)/(dN/dy) with T c for 
three values of g e ff- The thick solid horizontal lines rep- 
resents experimental value of (dN^, /dy) / (dN/dy) with 
associated error bars. 



ergies [53, l28j). The analysis of the experimental 
data reveals a lower bound on the transition temper- 
ature and the effective statistical degeneracy. The 
extracted value of T c ~ 162 ± 8 and the lower bound 
on g e f f ~ 18. The study of K/ir in the similar frame- 
work is under progress p9| . 
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